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^ ■ ABSTRACT 

^ . SMA observations of the massive star-forming region IRAS 18089-1732 in the 

1 mm and 850 /xm band reveal outflow and disk signatures in different molecular 

^ . lines. The SiO(5-4) data show a collimated outflow in the northern direction. In 

^ . contrast, the HCOOCH3 (20-19) line, which traces high-density gas, is confined 

i/^ , to the very center of the region and shows a velocity gradient across the core. 

^ ! The HCOOCH3 velocity gradient is not exactly perpendicular to the outflow 

Tj" ! axis but between an assumed disk plane and the outflow axis. We interpret these 

^ ! HCOOCH3 features as originating from a rotating disk that is influenced by the 

p ^. outflow and infall. Based on the (sub-)mm continuum emission, the mass of 

Q ! the central core is estimated to be around 38 Mq. The dynamical mass derived 

^ ! from the HCOOCH3 data is 22 Mq, of about the same order as the core mass. 

c§ ! Thus, the mass of the protostar/disk/envelope system is dominated by its disk 

j> ! and envelope. The two frequency continuum data of the core indicate a low dust 

^ , opacity index (3 ~ 1.2 in the outer part, decreasing to ~ 0.5 on shorter spatial 

^ ! scales. 
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submillimeter - techniques: interferometric 



-'^ Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138 

^Academia Sinica Institute of Astronomy and Astrophysics, No.l, Roosevelt Rd, Sec. 4, Taipei 106, 
Taiwan, R.O.C. 



- 2 - 



1. Introduction 

Unambiguous proof for disks in massive star formation is still missing. Millimeter contin- 
uum observations suggest flattened structures without providing velocity information (e.g., 
Shepherd et al. 2001), and molecular line studies suggest rotational motions but are often 
confused outflows and ambient gas (e.g., Zhang et al. 1998 and Beuther et al., this volume). 
Maser studies show disk signatures in some cases but are mostly not unambiguous as well 
(e.g., Churchwell 2002). The best evidence yet for genuine disk emission comes from CH3CN 
observations in IRAS 20126+4104 (Cesaroni et al. 1999). In this case, the velocity gradient 
defining the presence of the disk is aligned perpendicular to the bipolar outflow, consistent 
with the common disk/jet paradigm. To further investigate possible disk emission and its 
association with molecular jets, we used the Submillimeter Array (SMA) to observe the jet 
tracer SiO(5-4) and the hot-core tracer HCOOCH3 (20-19) in a massive star- forming region. 

The source IRAS 18089-1732 is a young High-Mass Protostellar Object (HMPO) which 
has been studied in detail over recent years. The source is part of a sample of 69 HM- 
POs selected mainly via infrared color-color criteria and the absence of strong cm emission 
(Sridharan et al. 2002). IRAS 18089-1732 is approximately at a distance of 3.6 kpc^ and its 
bolometric luminosity is about 10^-^ L© (Sridharan et al. 2002). Millimeter continuum obser- 
vations reveal a massive core > 2000 M© with H2O and CH3OH maser emission, and a weak 
1 mJy source is detected at 3.6 cm (Beuther et al. 2002a,c). As part of a single-dish CO out- 
flow study, wing emission indicative of molecular outflows was detected but the CO map was 
too confused to define a bipolar outfiow (Beuther et al. 2002b). During these observations, 
Beuther et al. (2002b) also observed SiO(2-l) at 3 mm, and bipolar structure was detected 
in the north-south direction. Furthermore, Sridharan et al. (2002) reported the detection 
of the hot-core-tracing molecules CH3CN and CH3OH. This letter focuses on the jet/disk 
observations and the (sub-)mm continuum data. A description of the line forest observed 
simultaneously is presented in an accompanying paper (Beuther et al., this volume). 



^The kinematic distance ambiguity is solved by associating the region via the near- and mid- infrared 
surveys 2MASS and MSX on larger scales with sources of known distance (Bontemps, priv. comm.). 
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2. The Submillimeter Array (SMA) 

IRAS 18089-1732 was observed with the SMA"^ between May and July 2003 in two 
different configurations with 3 to 5 antennas in the array. The phase reference center of the 
observations was R.A.[J2000] 18:11:51.4 and Dec.[J2000] -17:31:28.5. The frequency was 
centered on the SiO(5 4) fine at 217.105 GHz, the HCOOCH3 (20-19) line at 216.967GHz 
could be observed simultaneously in the same band. The HCOOCH3 fine consists of 8 distinct 
components but is dominated by 4 of them which are separated by 2.5 MHz (corresponding 
to 3.5 km s^^). The correlator bandwidth at that time was 1 GHz with a frequency resohition 
of 0.825 MHz. We smoothed the SiO(5-4) data to a spectral resolution of 3kms~^ and the 
HCOOCH3(20-19) data to 2kms~^ to increase the signal-to-noise ratio. The continuum was 
constructed via averaging the line-free channels in the upper side-band. The beam size at 
217 GHz was 2.7" x 1.7" and at 354 GHz 1.4" x 0.9". System temperatures in the 850 /xm 
band were between 300-900 K and in the 1mm band around 200 K. The continuum rms at 
217 GHz was ~ 8mJy and at 354 GHz 40mJy. The flux cahbration was estimated to be 
accurate to 25%. For more details on the SMA, the observations and data reduction, see the 
accompanying papers by Ho, Moran & Lo and Beuther et al. (this volume). 



3. Results 

3.1. Dust continuum emission 

Figure 1 compares the (sub-)mm continuum observations and shows additional cm con- 
tinuum and H2O and CH3OH maser data (Beuther et al. 2002c). Even in the highest- 
spatial-resolution data at 850 /xm, the dust emission remains singly peaked, i.e., it does 
not split up into multiple sources as observed in other massive star-forming regions, e.g., 
IRAS 19410+2336 (Beuther et al. 2003). Nevertheless, in our 1mm data we resolve elon- 
gated emission in the south and north-west, which demonstrates that IRAS 18089-1732 has 
a compact mm core with extended halo emission (Fig. 1). The halo emission is not seen in 
the 850 //m observations because of the reduced sensitivity and uv-coverage. While the weak 
3.6 cm peak and the H2O maser position coincide exactly with the (sub-)mm continuum 
peak, the CH3OH maser position is about 1.4" to the south. The latter could indicate that 
there might be a second source at the position of the CH3OH maser which we cannot dis- 



"^The Submillimeter Array is a joint project between the Smithsonian Astrophysical Observatory and the 
Academia Sinica Institute of Astronomy and Astrophysics, and is funded by the Smithsonian Institution and 
the Academia Sinica. 
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tinguish. Table 1 shows the derived peak and integrated fluxes (-S'peak and Sint) at 1 mm and 
850 /xm. Comparing the SMA 1 mm data with single-dish observations of the region (Beuther 
et al. 2002a), we find that about 85% of the fiux is filtered out in the interferometric data. 

It is difficult to derive a spectral index from the continuum images because the different 
uv-coverages filter out different amounts of fiux. However, we can measure fiuxes Suy in the 
uv-plane. Ideally, one would select the same regions in the uv-plane, but as this would 
reduce the amount of available data even more, it is reasonable to compare the values 
for matching baseline ranges (in units of k\). We selected one range of short baselines 
(20 — 40 k\, corresponding to spatial scales between 10.7" and 5.4") and one range of longer 
baselines (60 — 80 kX, corresponding to spatial scales between 3.6" and 2.7") where there were 
sufficient data in both frequency bands: the fiux values are shown in Table 1. The 3.6 cm 
fiux is only 0.9 mJy (Sridharan et al. 2002), and assuming free- free emission its contribution 
to the sub-(mm) observations is negligible. Assuming a power-law relation S oc u^'^^ in 
the Ray leigh- Jeans limit with the dust opacity index j3, we find f3 ~ 1.2 for short baselines 
corresponding to large spatial scales and (3 ~ 0.5 for large baselines corresponding to small 
spatial scales. These values are lower than the canonical value of 2 (Hildebrand 1983) 
and decrease to small spatial scales. The exact values should be taken with caution due 
to the large calibration uncertainty (estimated to be within 25%), but the trend of a low 
(3 decreasing with decreasing spatial scales appears real. Without mapping the selected 
sets of uv-data in the image-plane, we cannot determine whether the smaller spatial scales 
correspond to the central core or to another unresolved structure. However, the 850 //m 
continuum image, which is based on larger baselines, is more compact than the 1.3 mm 
emission (Fig. 1) indicating that the largest baselines do trace the emission from the core 
center. A low value of f3 has also recently been found in another high-mass star-forming 
region (Kumar et al. 2003), and Hogcrhcijdc & Sandell (2000) report that f3 decreases in 
L1489 from 1.5-2 in the envelope down to i3 = 1 at the inner peak. This trend may be due 
to grain growth within the central core/disk or high optical depth. 

Assuming optically thin dust emission at mm wavelength, we calculate the mass and 
peak column density using the 1.3 mm data following the procedure outlined for the single- 
dish dust continuum data by Beuther et al. (2002a). We use a grain emissivity index (3 = 1 
and a temperature typical for hot-core-like sources of 100 K. The central core mass then is 
~ 38 M© and the peak column density ~ 8 x 10^^ cm~^, corresponding to a visual extinction 
74v[mag] = NH2/(10^^cm'~^) ~ 800. As discussed by Beuther et al. (2002a), the errors are 
dominated by systematics, e.g., exact knowledge of (3 or the temperature. We estimate the 
masses and column densities to be correct within a factor of 5. 
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3.2. The outflow in SiO 

Figure 2 presents a channel map of SiO (5-4), and we find blue and red emission (systemic 
velocity 34.9 km s~^, §3.3) north of the continuum core. The SiO emission by itself allows 
an interpretation of a bipolar outflow centered 2" — 3" north of the core. However, we do 
not detect any sub-mm continuum source there (850 /im la mass sensitivity ~ 1 Mq), and a 
massive outflow without a sub-mm continuum source at the driving center is very unlikely. 
Therefore, we favor the interpretation of an outflow emanating from the core toward the 
north with a position-angle (P.A.) of ~ 20°. The fact that we see blue and red SiO emission 
toward the northern lobe indicates that the outflow axis is near the plane of the sky. For the 
red lobe, we find an increase in velocity with distance from the center of the core resembling 
the Hubble-law within molecular outflows (e.g., Downes & Ray 1999). 

Previous unpublished single-dish SiO (2-1) observations with 29" resolution also show a 
bipolar outflow in the north-south direction. The main difference is that the southern part 
of the SiO(2-l) outffow is not observed at the higher frequency and higher spatial resolution 
with the SMA. We can attribute this difference to several possible reasons. First, interfer- 
ometers fllter out the large scale emission, and without additional short spacing information, 
only the more compact jet-like emission is detected with the SMA. Second, the excitation 
of the 5-4 transition is higher with respect to the 2-1 transition. While the SiO(2-l) fine 
traces the lower temperature gas, the SiO (5-4) line is sensitive to the more excited gas being 
associated with the more coUimated component of the outflow. Furthermore, an asymmetry 
in the small-scale distribution of the dense gas can also contribute to the differences of the 
large-scale SiO (2-1) and small-scale SiO (5-4) observations. 

3.3. HCOOCH3 tracing the disk? 

The high-density-gas tracing molecular HCOOCH3 (20-19) emission is conflned to the 
central core traced by the (sub-)mm continuum (Fig. 2). A Gaussian flt to the central 
spectrum results in a systemic velocity of ~ 34.9 km s~^, and we observe a velocity gradient 
across the core. It is possible to fit the peak position of each spectral channel shown in Figure 
2 to a higher accuracy than the nominal spatial resolution, down to 0.5 HPBW/(S/N), with 
S/N the signal-to-noise ratio (Reid et al. 1988). We performed these fits in the uv-plane 
to avoid artifacts due to inverting and cleaning the data, the resulting positions are shown 
in Figure 3. The derived velocity gradient follows a P.A. of ~ 55°(±10°) from blue to red 
velocities. The separation of the most blue and most red positions is ~ 0.6", corresponding 
to ~ 2200 AU. Figure 3 also presents a position-velocity diagram for the fitted channels at 
the P.A. of 55° conflrming the space-velocity shift. This velocity gradient is neither in the 
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direction of the SiO outflow nor directly perpendicular to it. 

A velocity gradient solely due to a rotating disk would be perpendicular to the outflow. 
The SiO outflow is near the plane of the sky complicating an accurate determination of its 
P. A., and the outflow might also be affected by precession (e.g., IRAS 20126+4104, Shepherd 
et al. 2000). Regarding the HCOOCH3 emission, it is unhkely that this high-density tracer 
depicts a different outflow which would remain undetected by the SiO. Therefore, we regard 
neither the partially uncertain outfiow P.A. nor a second outflow sufficient to explain the 
difference in P.A. between the SiO and the HCOOCH3 emission. However, as shown in the 
accompanying paper (Beuthcr ct al., this volume), most molecular line data arc influenced 
by the molecular outflow, and it is likely that HCOOCH3 is also affected by the outflow. 
Furthermore, Ohashi et al. (1996) have shown that infall can influence the dense-gas velocity 
signature as well. Assuming that the inner core contains a rotating disk, which has an orien- 
tation perpendicular to the outflow, and that the HCOOCH3 velocity signature is influenced 
by the disk, the outflow and infall, it is plausible that the observed axis of the HCOOCH3 
velocity gradient is offset from the expected disk orientation. We do not observe the inner 
region of the disk where the outflow/jet is accelerated, but we may have detected the outer 
parts of the accretion-disk. 

Assuming equilibrium between centrifugal and gravitational forces at the outer ra- 
dius of the disk, we calculate the dynamical mass enclosed within the central 0.6" to be 
~ 22/sin^iMQ (with i the unknown inclination angle between the disk plane and the line of 
sight, and an HCOOCH3 velocity shift^ of 4.25 km s^^), of the order of the core mass derived 
from the dust emission (Table 1). Thus, in contrast to T Tauri systems where the disk mass 
is neghgible compared with the protostellar mass, in IRAS 18089-1732 a considerable fraction 
of the enclosed mass seems to be part of a large accretion disk and/or rotating envelope. 
This result flts the picture that IRAS 18089-1732 is in a very young evolutionary state, with 
the central (proto)star still accreting material from the surrounding core/disk. 

4. Discussion 

The combined SiO (5 4) and HCOOCH3(20-19) observations toward IRAS 18089-1732 
support a massive star formation scenario where high-mass stars form in a similar fashion as 
their low- mass counterparts, i.e., via disk accretion accompanied by coUimated jets/outflows. 
The HCOOCH3 observations still barely resolve the disk/envelope structure and the inter- 



^The full-width zero intensity is 12kms ^. Taking into account the 4 dominant HCOOCH3 line compo- 
nents spanning S.Skms'^ (§2), the HCOOCH3 velocity shift is (12 - 3.5)/2 = 4.25. 
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pretation is not entirely unambiguous, but the data indicate rotation which might stem at 
least partly from an accretion disk. Higher resolution observations in different disk tracers 
are needed to investigate the disk/envelope conditions in more detail. Furthermore, the con- 
tinuum data indicate a lower (3 than the standard value of 2, and we observe a decreasing 
j3 with decreasing spatial scales. As nearly all the gas observed in dust emission takes part 
in the assumed dynamical rotation observed in HCOOCH3 this low /3 might be due to grain 
growth or high opacity within a disk-like structure (e.g., Beckwith et al. 2000). 

It's a pleasure to thank the referee Riccardo Ccsaroni for his comprehensive comments 
improving the quality of the paper. We also like to thank a lot the whole SMA staff for mak- 
ing this instrument possible! H.B. acknowledges financial support by the Emmy-Noether- 
Program of the Deutsche Forschungsgemeinschaft (DFG, grant BE2578/1). 
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Tablc 1. (STib-)millimctcr continmim data 
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'Uncertainties are discussed in the main text. 
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Fig. 1. — The grey-scale (with dotted contours) shows the SMA 1.3 mm continuum emission 
in both panels. The left panel shows in heavy contours the SMA 850 /xm continuum emission 
and the 850 /xm synthesized beam. In the right panel we also present the cm continuum 
emission and H2O and CH3OH maser position (taken from Beuther et al. 2002c), and the 
1.3 mm synthesized beam. 
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Fig. 2. — SiO(5-4) (top) and HCOOCH3 (bottom) channel maps. The grey-scale outlines 
the 850 nm continuum emission. At the bottom left of each panel, the filled ellipses shows 
the 850 /xm beam and the unfilled ellipses the 1.3 mm beam. LSR velocities (in kms~^) are 
given in the upper left corner of each panel. The LSR velocities for HCOOCH3 are calculated 
with respect to the line component at 216.9662 GHz, the systemic velocity is ~ 34.9 km s~^. 
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Fig. 3. — Presented are the position-position (p-p) and position-velocity (p-v) diagrams of 
HCOOCH3. The offsets in the p-p plot were fitted in the uv-plane for each velocity channel 
presented in Fig. 2. The positions delineate the P.A. of 55°, and the error bars represent 
statistical errors of the fits. The open pentagon marks the position of the 1.3 mm continuum 
peak. The offsets in p-v diagram are offsets from the continuum center position along the 
assumed disk plane (P.A. 55°). 



